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In this talk, we describe part of a recent work [T] on the correlation function of a, V — A current with slV + A 
current in the framework of QCD in the limit of a large number of colours Nc- The discussion takes place within 
two successive approximations of this theory, called MHA and MHA+V. Results concerning the evaluation of 
chiral condensates of dimension six and eight, as well as matrix elements of the Q7 and Qg electroweak penguin 
operators, are given. 



1. Introduction 

The correlation function of a left-handed cur- 
rent with a right-handed current 

= u(a;)7^i^d(x) 

and 

i?^(o) = dW^^^W, 

in QCD and in the chiral limit, depends only on 
one invariant amplitude IIlr{Q^) of the euclidean 
momentum squared = —q^, with q the mo- 
mentum flowing through the two-point function: 

2i yd*a;e*«-^(0|T{L^(x) -R^IO)} |0) 

= (q'-q'' ~ q'gn n^RiQ') . (1) 

Here, we shall be particularly concerned with 
the study of IlLii{Q^) in the limit of a large 
number of colors Nc in QCD. We want to com- 
pute the lowest dimension operator product ex- 
pansion (OPE) condensates of this function in 
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two successive approximations to the large-A''^ 
limit: the so-called minimal hadronic approxima- 
tion (MHA) P] consisting of a spectrum of the 
pion state, a vector state and an axial vector 
state; and the improved approximation where an 
extra higher vector state is added. 

Part of our motivation comes from the fact that 
in the literature the values of these condensates 
are rather controversial. Phenomenological anal- 
yses in refs. |2I3I4I5| find for example opposite 
sign for the first two condensates, in contradis- 
tinction to the results in refs. |6l7l8j which find 
the same sign. 

In this talk, we give results only for the two 
lowest dimension condensates and the matrix el- 
ements of and Q% (see however for more 
results and details). 

In practice, we shall be working with the di- 
mensionless complex function Wlr[z] defined as 



Wlb\z] = -zUlr{zMp) , with Rez = 




and use the mass of the lowest massive state, the 
p{770 MeV), to normalize quantities with dimen- 
sions. In large- A'^c QCD the function Wlr[z] is a 
mcromorphic function and, therefore, in full gen- 
erality, it can be approximated by successive par- 
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tial fractions of the type 

p ^ N 
Wlb.[z]^AnI[^—-1[{z + <Jj), (2) 

with pi — \ and pi 7^ pk for i ^ k\ where P 
(and A^) get larger and larger, but finite. Ajv is 
the overall normalization and in what follows, we 
shall often use the notation p2 = pa and p^ = 
Pv- 

On the other hand, in QCD, the OPE of 
the two currents in Eq. fixes the large-Q^ 
fall off in 1/ Q^-powers of the invariant function 

nLfl(Q2) [Hj to 

00 



n=l 



x(0.„,4A^^)>(Q^)-<"^^^-^f:py^. (3) 

Matching the leading asymptotic behaviour for 
large-z in Eq. (01 to the one of the OPE in 
Eq. restricts the number of zeros N and the 
number of poles P in Eq. ^ to obey the con- 
straint 



N - P 



(4) 



The case where iV = corresponds to the MHA 
while the MHA+V case is obtained for = 1. 

2. T\lr in the MHA and t:-V-A~V' Spec- 
tra 

It is claimed by some of the authors of refs. 
I7I8| that the reason why their phenomenological 
analysis of the chiral condensates give the same 
sign for (Og) and (Og) is due to the fact that 
the hadronic r-decay spectrum is sensitive to the 
presence of the p', while the MHA ignores all 
higher states beyond the first axial state. Partly 
motivated by this claim, we want to analyze here 
the case, beyond the MHA, where an extra vec- 
tor state y, and therefore one zero cr, are also 
included. Let us collect some relevant equations 
describing this case (the corresponding expres- 
sions in the MHA case can be obtained by putting 
Pv = o- in the MHA-I-V equations). 



2.1. The Correlation Function 

With a spectrum of the pion pole, and V, A, 
and V' states, the relevant correlation function is 



Wlr[z]^Ai 
where Ai 



z + a 



- PAPv 

-2 



[z + l){z + pa){z + pv) ' 
= Pf and therefore 
2 



— = 



a 



-My-pFPAPV ■ (5) 



Because of the positivity of Wlr[z] for Re 2; > 0, 
the position of the zero has to be in the Minkowski 
axis and, then, a > 0. 

2.2. The Linear Constraint 

This is a very interesting constraint, key equa- 
tion of our analysis, which already provides a 
semi-quantitative argument in favor of the op- 
posite sign option for the condensates {Oq) and 
(Os)- It simply follows by expanding Eq. Q to 
first non-trivial order in inverse powers of z 



N 



P 



1 m 

M2 (Oe) 



(6) 



In the case corresponding to the MHA, where by 
definition there are no zeros, this constraint sim- 
ply becomes: 



1 



PA 



1 (Os) 



(7) 



implying that, in the MHA (Og) and (Og) must 
have opposite signs. 

Now, in the MHA+V case, Eq. ((HJ reduces to 



a - {1 + PA + Pv) = 



1 (Os) 



(8) 



This equation states that for (Og) to have the 
same sign as (Og), the position of the zero has to 
be far beyond the largest V -pole: 



a > 1 + PA + Pv 



(9) 



Fixing the position of the poles at the values of 
the observed spectrum (and ignoring errors for 
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the purpose of the discussion), one has My = 
0.776 GeV, Ma = 1.230 GeV {pA = 2.5) and 
Myr = 1.465 GeV {py> = 3.6); which means that 
for the equal sign requirement option to be sat- 
isfied, one must have a > 7.2. In GeV units this 
corresponds to a mass of 2.1 GeV. Now, in writ- 
ing a large- A'c ansatz for the Wlr[z] function, 
one is implicitly assuming an effective cancella- 
tion between the extra poles and zeros in the com- 
plex z-plane which lie beyond a disc of radius sq 
covering all the poles and zeros retained in that 
approximation. The result cr > 7.2 implies that 
the radius in question has to be ^/sq > 2.1 GeV. 
A priori that seems a good thing because the 
OPE-matching is now applied at > sq; i.e. 
in a more asymptotic region than in the case of 
the MHA ansatz; however, it also implies that 
there are no further poles in the region between 
My I ~ 1.5 and the effective mass M^ ~ 2.1 GeV 
corresponding to the zero at tr ~ 7.2. This, how- 
ever, is in contradiction with the observed oi- 
like state at Ma> — 1.64 GeV and p-like states 
at Mv>, ~ 1.72 GeV and My,„ ~ 1.9 GeV be- 
low Afg- ~ 2.1 GeV. Alternatively, if one ex- 
cludes those three states A', V" and V" as all 
the phenomenological analyses using r~data do 
in fact, then the position of the zero a should be 

< -^r- ~ 4.5, implying according to Eq. 0, 
that (08 ) and (Oe) must have opposite signs, in 
contradiction with the claims of refs. |6I7I8 || . 



3. Numerical Analyses and Conclusions 

Confronting the MHA and MHA-I-V' approxi- 
mations to the experimental values of some ob- 
servables that we introduce in the following will 
allow us to test their consistency and adjust their 
free parameters {pp and pA in the MHA case, pp, 
PA, p'y and a in the MHA-I-V case) in order to 
make predictions for the condensates and matrix 
elements. We use as input the following set of 
experimental data 



•TTJ 



(640 ± 246) X 10"^ MeV, 



ref. [TT] 
ref. 

= (6.77 ± 0.32) X 10-3 MeV , ref. [TT] 



6m„ = 4.5936 ± 0.0005 MeV , 
Lio = (-5.13 ±0.19) X 10~3, 



Lg = (6.9 ±0.7) X 10"^ 
Mp = (775.9 ±0.5) MeV; 



ref. HH 
ref. [12 
ref. [TT] 



In fact, some of these observables, when ex- 
pressed in terms of the MHA or MHA±V' pa- 
rameters, depend not only on these ones but also 
on (m^++m^o) and/or My. Therefore, it is more 
appropriate for our purposes to use in our fit pro- 
cedure the dimensionless quantities: 



m^++m„o 3apA'^og{pA) 



Ml 



4 TT PA - 1 



Lro = ~ip.(l + ^), 
J_ _ 47ra^ PA 

Mp''^^'^^- 24PA-1' 



3 PA - 1 ' 
a y/pA 



(10) 

(11) 
(12) 

(13) 
(14) 



For lack of place we gave the formulae only in the 
case of the MHA, the corresponding MHA+V 
equations can be found in pp. 

A standard statistical regression method 
permits to fit the parameters. 

3.1. The case of the MHA spectrum 

For this spectrum, we find pp = (12.36±0.35) x 
10-3 and PA = l-464±0.004,withaxmin = 1-21 
for 3 degrees of freedom (dof). The covariance 
matrix is given by 



, . /1.21 

cov(pf,pa) = L 3g 



1.36^ 

162 ; 



X 10"^ . (15) 



We deduce from these results two conclusions: 
first that the MHA framework is statistically rel- 
evant and second that the fitted free parameters 
have small statistical errors. 

3.2. The case oian— V~A — V' spectrum 
Eq. H14|l becomes now a function of gy as ex- 
plained in The way we treat this is by con- 
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sidering the observable^: 



1 

1 a{pF - 



2Lgpv') (pa - - 1) 



(1 - pv'YpfPaPv (o- - 1) 



as a function of Lg which has an error itself, and 
is added as an extra parameter in our fit. The 
number of dof does not change since Lg is also 
taken as an observable. For this case, we impose 
a criterion of rejection through the ordering: 

So 



PA < pv < < Po 



The first and second inequalities reflect the 
knowledge that the new state has a higher mass 
than the axial, that it is a y-like pole and, 
therefore, its residue contributes positively to 
the WL/{(2:)-function; the third inequahty fol- 
lows from the requirement that the perturbative 
threshold Sg already lies beyond the radius where 
the analytic structure of the poles and zeros re- 
tained satisfies the leading OPE constraint. 

The regression leads to 



PF = (12.36 ±0.03) X 10" 
PA = 1.466 ± 0.003, 
pv = 2.63 ±0.01, 
a = 2.64 ±0.01. 



(16) 
(17) 
(18) 
(19) 



The results in Eqs. 1)16(1 to p9(l correspond to a 
value: Lg = (6.44 ± 0.02) x lO^^, with a xLn = 
0.60 for 1 dof. 

We conclude that the parameters pp and pA are 
statistically stable when compared to those found 
in the MHA case. We also find that pv' ~ c, 
which is consistent with the fact that the MHA 
approximation already seems to have the bulk of 
the full large-TVc information. In other words, 
adding an extra V'-pole appears to be compen- 
sated, at a very good approximation, by the po- 
sition of the nearby zero. 

3.3. predictions 

Using Eqs. (jSJ, lO (and corresponding equa- 
tions in the MHA case) as well as expressions 

^The decay width of p ^ tttt is Tp^^^ = (150.4±1.3)MeV 

ini. 



of the matrix elements of four-quark operators 
mM), we find 





xlO^ GeV^ 


xlO^ GeV* 


MHA + V 


-7.90 ±0.20 ±1.62 


+11.69 ±0.32 ±2.53 


MHA 


-7.89 ±0.23 ±2.01 


+11.71 ±0.34 ±3.08 





M7 (NDR) 


A'h (HV) 


MHA + V 


0.12 ±0.00 ±0.01 


0.59 ±0.01 ±0.06 


MHA 


0.12 ±0.00 ±0.02 


0.59 ±0.01 ±0.11 




Ms (NDR) 


Afg (HV) 


MHA + V 


2.00 ±0.03 ±0.20 


2.15 ±0.03 ±0.22 


MHA 


1.99 ±0.03 ±0.36 


2.15 ±0.03 ±0.39 



The way of obtention of these results as well 
as the origin of systematic errors are explained in 

n 

We then see that within errors, the two sets of 
predictions from MHA and from MHA + V are 
perfectly consistent with each other. 
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